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Abstract 

The investigation of a sizable thermal enhancement of magnetization is put forward for uniaxial 
ferromagnetic nanoparticles that are placed in a rotating magnetic field. We elucidate the nature 
of this phenomenon and evaluate the resonant frequency dependence of the induced magnetization. 
Moreover, we reveal the role of magnetic dipolar interactions, point out potential applications and 
reason the feasibility of an experimental observation of this effect. 
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Presently, the study of magnetic nanoparticles and their structures is one of the most 
important research areas in nanoscale physics. A first reason is that such nanoparticles in- 
creasingly find numerous applications that range from medicine to nanotechnology. Another 
reason is that these systems exhibit a number of remarkable physical phenomena, such as 
quantum tunneling of magnetization giant magnetoresistance I^J, exchange bias and 
finite-size and surface effects Q], to name but a few. Moreover, the study of fundamentals 
of magnetic behavior in these systems is also an important issue, especially for high-density 
data storage devices 

From a practical point of view, the lifetime of stored data and the switching time (i.e., 
the time during which the reversal of the nanoparticle magnetic moments occurs) are salient 
characteristics of such devices. Now, thanks to the experimental discovery of fast switching 
of magnetization |6'] , the switching time reaches the fundamental (picosecond) limit for field- 
induced magnetization reversal. On the contrary, a feasible lifetime must cover up to ten 
years and beyond. Its value is usually limited by the superparamagnetic effect jlj and is 
defined by the probabilities that the nanoparticle magnetic moment m stays in the up 
(cr = +1) and down (a = —1) equilibrium directions. These probabilities, which are also 
responsible for other thermal effects in such systems including magnetic relaxation Q], are 
very sensitive to small perturbations that change the static states of the magnetic moments. 
Namely, according to the Arrhenius law |9] the ratio p+i/p-i is approximately given by 
exp{AE / kT) , where AE = Ej^i — E^i, E^r is the potential barrier for the reorientation a — *■ 
—a, k is the Boltzmann constant, and T is the absolute temperature. Therefore, if without 
perturbations AE = 0, then = 1 and the nanoparticle system is demagnetized. 

But due to the exponential dependence on AE and T, the ratio p+i/p-i can drastically be 
changed by small perturbations. In particular, a static magnetic field H applied along the 
nanoparticle easy axis of magnetization yields AE = 2Hm (m = |m|), and so p+i strongly 
differs from p_i if \H\/Ha ^ l/4a where a = Ham/2kT, and Ha is the anisotropy field. 
This means that even small magnetic fields (in comparison with Ha) almost fully magnetize 
the nanoparticle systems when a ^ 1. 

In the case of time-periodic perturbations the situation is not settled yet and far less re- 
searched. On the one hand, these perturbations generate dynamical states of the nanoparti- 
cle magnetic moments that, because of their natural precession, are expected to be different 
for the up and down magnetic moments. Therefore, a dynamical magnetization of the sys- 
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tern, i.e., a magnetization which is induced by periodic perturbations at T = 0, can exist but 
it is expected to be small and, at first sight, it cannot be changed by thermal fluctuations 
since these perturbations do not change the mentioned above potential barriers. On the 
other hand, it may be expected that, due to the different precessional states of the up and 
down magnetic moments, the probabilities p^j are different and thus thermal fluctuations 
contribute to the dynamical magnetization. In this Letter, we attempt to solve this chal- 
lenge which comprises both basic and applied aspects of the concept of a mean first-passage 
time. 

Dynamical magnetization. — To calculate the dynamical magnetization, we consider the 
case of identical, non-interacting nanoparticles whose easy axes of magnetization are parallel 
to each other and the dynamics of the magnetic moment m is governed by the deterministic 
Landau-Lifshitz equation jlo| 

rri = —7m X He// m X (m x He//). (1) 

Here 7(> 0) is the gyromagnetic ratio, A(> 0) is the dimensionless damping parameter. 
He// = —dW/dm is the effective magnetic field acting on m, and W is the nanoparticle 
magnetic energy. If the easy exes are parallel to the z axis and the external magnetic 
field h{t) is circularly polarized in the xy plane, i.e., h{t) = h{cos ujt, p sin out, 0), where 
h = \h{t)\, uj is the frequency of field rotation, and p = — 1 or +1 (the sign — corresponds 
to the clockwise rotation of h(t) and the sign + to the counterclockwise one), then 

W = \mHa sin^ 6 — mh sin 6 cos ip (2) 

[ip = if — put, and 6 and ip are the polar and azimuthal angles of m, respectively). For this 



case, the solution of Eq. (HI) is well studied in the context of ferromagnetic resonance 111 

n 

and nonlinear magnetization dynamics Specifically, the precession angles 6^ of the up 
and down magnetic moments (see Fig. 1) at ^ 1 are given by 

^ , ^ (3) 

where uJr = '~fHa. Note that the angles 6cr exhibit the resonance dependence on u only for 
pa = +1. 

We define the dimensionless magnetization of the nanoparticle system induced by the 
magnetic field h(t) as p = {l/N) Y^j^i rrizi/m (the index i labels the nanoparticles, ^ 1). 



If the magnetic anisotropy barrier essentially exceeds the thermal energy, i.e., a S> 1, this 
definition yields fi = {l/N)'}^^aN(jCos9fj [No- is the number of magnetic moments in the 
state cr, A^_i + A^+i = A^), and so /i = X^o- "^Pcr cos 6*0- ijPa = ^a/N). Assuming = 1/2 at 
T = (this means that the nanoparticle system is demagnetized if /i = 0) and using the 
condition 9a- 1 and Eq. (jS)), for the dynamical magnetization fi^ = ^\t=o = {d'^i — ^+i)/4 
we obtain 

~ ^[(l + A2V2+a;2]2_4^2^2- i ) 

According to this result, the direction of dynamical magnetization and the direction 
of magnetic field rotation follow the left-hand rule (the reason is that the natural precession 
of the magnetic moments is counterclockwise), and the dependence of on uj always has a 
resonant character with max = ^d\u)=uimi where ujm = (co'r/-\/3)[l — A^ + 2(l + A^ + A'')^/^]"'^/^. 
But the value of ^d\uj=LOm is very small, however, because with max^^o- = h/XHa even for 
A ^ 1, Eq. (jH) yields ^d\uj=LUm = —p{h/2XHaY ■ Note also that after switching on the 
magnetic field h(t) the initially demagnetized system reaches the steady-state magnetization 
yUrf during a time interval of the order of tr = 2/XuJr (we recall that for T = the states a 
of the magnetic moments are not changed with time and = 1/2). 

Thermal enhancement of the dynamical magnetization. — If T 7^ then the dynamics of 
the magnetic moments becomes stochastic. In this case, due to thermal fiuctuations, the 
magnetic moments can perform random transitions from the one state a to the other —a 
and the probabilities Pa can thus depend on h(t). But, in contrast to a static magnetic 
field, a rotating field has no preferential direction and so it does not impact p^ directly. 
Nevertheless, the probabilities p+i and p_i must be different in the presence of h(t). The 
reason is that if the mean times t^ which the magnetic moments reside in the states a are 
much larger than the precession time 27r/ci;, then the up and down magnetic moments spend 
almost all time near the conic surfaces with the cone angles 6^1 and 6'_i, respectively. Since 
these angles are different, see Eq. Q, the times t^ must be different as well. Accordingly, 
because in the steady state Pa = to-/(t+i + t-i), we conclude that the probabilities Pa are 
also different and so thermal fluctuations in fact do contribute to the induced magnetization 

Using the conditions 9^ <^ 1 and a S> 1, from the definition of /i we obtain fi = fit + fJ-d 
(we neglect the term UtiS'^-i + ^+i)/4), where Ht = p+i — p~i is the desired contribution 
arising from the joint action of thermal fiuctuations and rotating field. According to the 
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above argumentation, the condition p^r < p^a holds if 6cr > O-a- This imphes that the 
thermal contribution /i^ always enhances the deterministic part /i^. Moreover, one expects 
the enhancement increases with decreasing temperature. It is important to note in this 
context that /i denotes the equilibrium magnetization which is established during a time 
interval of the order of the transition time ttr between the states a and —a (t^r ~ max ^ tr 
and ttr — i> C)0 as T — > 0). At these times the probabilities p+i and p-i are generally different 
and thus limr^o ^^ 7^ l^d- 

For determining the mean residence times to- which define the probabilities p^ and the 
magnetization /i^, we used the mean first-passage 

to our situation is well-founded because in the case of small rotating field the magnetic 
moments that are in the state a reach any point of the separatrix, which separates the up and 



down states, with almost the same probabi 
of the magnetic moments is Markovian 



ity density. Given that the stochastic dynamics 
the standard mean first-passage procedure is 



employed in order to account for the influence of the rotating field. Specifically, starting out 
from the two-dimensional backward Fokker-Planck equation ^ in the rotating frame, 
we succeeded to derive a mathematically tractable and physically transparent expression for 
the mean residence times: 

to = to exp[a(-^^ + a2Hes/Ha)], (5) 



where to = tj-y ir/Aaexp a is the mean time which the magnetic moment spends in the up 
or down state at h = 0, and ifefr = —ppjh'^/uj (|-f^cfr| ^ Ha, p ~ 1). 

As follows from (0), the rotating magnetic field infiuences the mean residence times via 
two different mechanisms. A first one consists in the appearance of the dynamical states of 
the magnetic moments which are characterized by the precession angles 6^- The contribution 
of these states to to- is governed by the first term in the argument of the exponential function, 
which always decreases to-. The second mechanism consists in changing the effective potential 
barrier between the up and the down states. Its contribution to to is described by the second 
term in the argument of the exponential function, where H^s can be interpreted as a static 
effective magnetic field applied along the easy axis of magnetization. Since the sign of Hes 
depends on p, this mechanism can either increase or decrease to. 

Using Eqs. © and (jHl), the definition fit = (t+i — t_i)/(t+i +t_i) leads to our main result 

fit = tanh[2a(/irf + H^s/Ha)]. (6) 



It shows that the rotating field always magnetizes the nanoparticle system perpendicular 
to the plane of field rotation and the direction of magnetization is uniquely defined by the 
direction of field rotation: sgn fi^o = —p- In general, both mentioned mechanisms contribute 
to fif However, in the most interesting resonant case, when uj ~ uj^ and A <^ 1, the first 
mechanism is dominating {fij^Ha/Hcs ~ A^^ ^ 1), and thus the magnetization (jH)) becomes 



According to this relation, being valid for a ^ 1 and \fid\ ^ 1, the condition pt/pd ^ 1 
always holds (specifically, if a\fid\ ^ 1 then fit/Pd ~ 2a). This means that a small dynamical 
magnetization is strongly enhanced by thermal fiuctuations, i.e., /i ~ fif Comparing (|7j) 
with the magnetization of an Ising paramagnet, tanh(miJ/A;T), we see that the magnetic 
field rotating in a plane perpendicular to the easy axes of the nanoparticles acts as a static 
magnetic field H = Hafid, which is applied along these axes. As in the case with fid, the 
induced magnetization fit as a function of u exhibits a resonance character. The dependence 
of fit on the reduced frequency uj = uj/uj^ is depicted in Fig. 2, curve 1, for a system of 
spherical nanoparticles with Ha = 6400 Oe, m/V = 1400 G {V is the nanoparticle volume), 
r = 4nm (r is the nanoparticle radius), A = 10~^, h = 10 Oe, p = —1, and T = 300 K. Note 
that for these parameters a = Ha'm/2kT ^ 29, pt\cj=i ~ 0.34, pd\d}=i ~ 6.1 x 10~^, and 
Pt/Pd ~ 56. 

Thus, the above results show that the magnetic field rotating in the plane perpendic- 
ular to the easy axes of magnetic nanoparticles changes the probabilities of the up and 
down orientations of the nanoparticle magnetic moments. Due to thermal fiuctuations, 
the nanoparticle system magnetizes in the direction with larger probability. In the case of 
nanoparticles with large anisotropy barrier (when a ^ 1) this contribution, i.e. pt, to the 
total magnetization p considerably exceeds a small dynamical contribution, i.e. pd, which 
is induced by the rotating magnetic field at T = 0. We emphasize that the magnetization 
of nanoparticle systems in the rotating magnetic field arises from the different dynamical 
behavior of the magnetic moments in the up and down states. In turn, the difference in the 
dynamics of the magnetic moments results from the existence of a well-defined direction of 
their natural precession (counterclockwise when viewed from above). 

Role of dipolar interactions. — To check the role of the magnetic dipolar interaction, 
we performed a Monte Carlo simulation for two-dimensional arrays of dipolar interacting 
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nanoparticles, representing an important class of patterned magnetic recording media j^. 
In doing so we assumed that the centers of nanoparticles occupy the sites of a square 
lattice of size Ld x Ld {L is a natural number, {L + 1)"^ = N, d is the lattice spacing). The 
easy axes of the nanoparticles are perpendicular to the lattice plane and the magnetic field 
rotates in this plane. In contrast to the previous case, the dipolar magnetic field acts on each 
magnetic moment. This field is changed from site to site and, due to the random motion of 
the magnetic moments, fluctuates with time. For a ^ 1, the fluctuations of the magnetic 
moments and the rates of their reorientations are small. Thus, the dipolar field acting on 
the i-th magnetic moment during the l-th step can be approximated as Hj(/) = (0, 0, Hi{l)), 
where Hi{l) = —mY,j^iO'j{l)/rl-, <Jj{l) = +1 or —1, and is the distance between the 
centers of the nanoparticles. In this case, the magnetization of the nanoparticle system can 
be represented through the step-dependent magnetization = Y.f=i o'i{l)/N as follows: 



{f^t) = j ,\, E/i.(0- (8) 



Here, to be sure that the system reached the steady state and the averaging procedure is 
correct, the conditions Zi ^ 1 and ^2 ~ ^ 1 are implied. 

In order to apply the Monte Carlo method for calculating (fit), we need to evaluate for 
all i and / the probabilities pcr.{i,l) that the ith magnetic moment stays in the states o"i(/). 
But in our case the conventional approach involving the Boltzmann factor for the solution 
of this problem is not applicable because the rotating field depends on time. Therefore, we 
extended the above method to the case of dipolar interacting nanoparticles and calculated 
the mean times that the ith magnetic moment spends in the up and down states: 

_ exp [a(cos6'^, + aibif - a] 
''^'"'^ (l-62)(cos^.,+aA) ' 

where bi = Hi{l)/Ha and 6^- is the precession angle of the ith magnetic moment that is 
defined by Eq. Q in which cr should be replaced by crj(/) and Ur by uJr + '~fHi{l) (for brevity, 
we omitted the arguments i and /). Next, defining PaXh = t„.{i, /) +t_i{i, I)] and 

using the numerical procedure developed in 3| , we performed a Monte Carlo simulation of 
the magnetization (fit) induced by the rotating field in two-dimensional systems of dipolar 
interacting nanoparticles. 

In Fig. 2, curve 2, we depict the dependence of {fit) on tD for the square array of the 
same nanoparticles driven by the same rotating field at iV = 10^, d = 5r, r] = 10~^, 



h — 10^, and Z2 = 5 x 10^. Comparing the curves 1 and 2 shows that the magnetic dipolar 
interaction reduces the induced magnetization, widens its frequency dependence, and raises 
the resonance frequency > !)• We emphasize that although the dipolar interaction 
reduces the induced magnetization, its experimental observation is still possible even in this 
strongly interacting case. Note also that fi^-^ , which we evaluated within the mean-field 
approximation (see Fig. 2, curve 3), distinctly differs from {/it) because this approximation 
does not account for the crucial feature of the dipolar interaction in these systems, i.e., its 
antiferromagnetic character. 

Potential applications. — The above results evidence that the frequency dependence of 
the induced magnetization is detectable and because the magnetic resonance methods are 
both very accurate and sensitive, its experimental determination thus provides valuable 
information about the dipolar field distribution in such systems. In particular, the average 
dipolar field acting on the resonant particles can approximately be estimated as pHa{Cjm~^)- 
Moreover, due to the selective change of the thermal stability of the magnetic moments, 
which is controlled by the characteristics of the rotating magnetic field, it calls for potential 
applications in magnetic recording technology. 

Resume. — Wc succeeded to show that a small dynamical magnetization of nanoparticle 
systems that is induced by the circularly polarized magnetic field is strongly enhanced by 
thermal fluctuations. The thermally enhanced magnetization exhibits as a function of the 
field frequency a resonant character, possessing a well pronounced extremum. The magnetic 
dipolar interaction increases the relative width of the frequency dependence of magnetization 
while causing a decrease of its strength and a corresponding shift of its maximum to higher 
frequencies. 
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FIG. 1: Sketch of the precession of the up and down magnetic moments (arrows show the directions 
of their natural precession) at p = — 1. 

FIG. 2: Plots of fit (curve 1), (/xj) (curve 2), and /i^^ (curve 3) as the functions of the reduced 
frequency a). 
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